Prospects of medium tomography using back-to-back hadron correlations 
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We discuss the prospects of extracting information about the bulk QCD matter distribution 
and evolution on the basis of hard hadronic back-to-back correlations in ultrarelativistic heavy- 
ion collisions. Using both hydrodynamical and parametrized evolution models for the spacetime 
evolution of the produced matter, which have been tested against RHIC data, we study six different 
setups for the spacetime dependence of hard-parton energy losses. Assuming that the energy loss of 
hard partons traversing the medium is radiative and calculable in the BDMPS formalism, we adjust 
one parameter, the quenching power scale, to the measured Raa in each of the setups and study 
the systematic variations of the back-to-back yield as a function of pr- We show which spacetime 
regions are probed by one-particle and two-particle observables and study the role of longitudinal and 
transverse expansion in some detail. We also comment on the importance of considering fluctuations 
around the average energy loss. We conclude that while current data are too limited in momentum 
coverage, future data for higher trigger energy might provide the lever arm in away side hadron 
momentum necessary to perform medium tomography, provided that sufficient precision can be 
achieved. 
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I. INTRODUCTION 



Announcements have been made by all four detector col- 
laborations at RHIC [l] that a new state of matter, dis- 
tinct from ordinary hadronic matter, has been created in 
ultrarelativistic heavy-ion collisions (URHIC). A new and 
exciting challenge for both experiment and theory is now 
to study its properties. The energy loss of hard partons 
created in the first moments of the collision has long been 
regarded a promising tool for this purpose [1,[1,I1,[1,[1,0|- 
So far, most of the effort in parton energy loss studies 
has been directed to understanding the nuclear suppres- 
sion factor Raa, i-e. the observed transverse momentum 
spectrum of hard hadrons divided by the scaled expecta- 
tion from proton-proton collisions. However, as recently 
argued in [8|], Raa exhibits only very limited sensitivity 
to the energy loss mechanism or properties of the medium 
beyond the fact that the quenching of jets is substantial. 
Thus, in order to overcome this obstacle and gain in- 
formation about the medium properties, one possibility 
outlined in (§] is to measure the momentum spectrum 
of hadrons correlated with a hard real photon at fixed 
energy. This has been shown to be able to discriminate 
between different scenarios of energy loss. 

In principle more detailed information about the en- 
ergy loss mechanism and the medium is also available in 
back-to-back hadron correlations. Measurements of two- 
particle correlations involving one hard trigger particle 
and associate hadrons above 1 GeV have shown a sur- 
prising splitting of the away side peak for all centralities 
but peripheral collisions, qualitatively very different from 
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a broadened away side peak observed in p-p or d-Au colli- 
sions [9( . Interpretations in terms of energy lost to prop- 
agating colourless [l(J, [ll|, [r| and coloured [l3[ sound 
modes have been suggested for this phenomenon, and cal- 
culations within a dynamical model evolution have shown 
that the data can be reproduced under the assumption 
that a substantial amount of lost energy excites a sonic 
Shockwave [lfj. Thus, it appears that for semi- hard as- 
sociate hadron momentum scales the recoil of the soft 
bulk medium is probed rather than the energy loss of 
the hard parton. While properties of the medium can 
be inferred from these measurements as well [l4j |. this is 
outside the scope of the present paper where we focus 
on the measurement of energy loss in hard back-to-back 
correlations. 

For sufficiently high associate hadron transverse momen- 
tum Pt > 4 GeV, back-to-back correlations with vacuum 
width are observed experimentally [TH, [l6| and the mea- 
sured yield per trigger is in agreement with the expec- 
tation from radiative energy loss in a dynamic medium 
[ItJ ■ As argued in Q , a measurement of Raa , or more 
general the suppression of a single hadron observable 
probes an averaged energy loss probability distribution 
(P(AE, E))t ab where the averaging is done over all pos- 
sible initial vertex positions (determined by the nuclear 
overlap Tab) and the soft matter. However, in two parti- 
cle correlations, the requirement that a trigger hadron is 
observed leads to a geometrical bias, and thus the yield 
per trigger of away side hadrons is determined by a differ- 
ent averaged energy loss probability (P(AE, E))t v where 
the averaging is done over all vertices leading to a trig- 
gered event. Thus, e.g. the distribution of away side 
pathlengths will be very different from the distribution 
underlying Raa ■ The question we would hence like to ad- 
dress is whether the difference between (P(AE, E))t ab 
and (P(AE, E))t t is (dependent on the medium evolu- 
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tion model) significant enough to infer non-trivial infor- 
mation about the medium, or in other words, to what 
extent there is information in back-to-back correlations 
beyond the information carried by Raa- 
Our strategy to investigate the question is as follows: 

1 . We consider two types of models for the medium evolu- 
tion which reproduce the observed bulk characteristics of 
Au-Au collisions at full RHIC energy ^snn = 200 GeV. 
Type I is the longitudinally boost-invariant hydrody- 
namical evolution model discussed in [lSj ]. where the 
multiplicities and transverse momentum spectra of pi- 
ons, kaons and protons in Au-Au collisions at ^/s^n — 
130 GeV are used as constraints and the initial condi- 
tions are computed from perturbative QCD+saturation 

19 1 . Type II is the parametrized evolution model of Ref. 

20] , where in addition to the Pt spectra also the rapidity 
distribution [2l[ and Hanbury-Brown-Twiss (HBT) cor- 
relations [22j] of central 200 AGeV Au-Au collisions at 
RHIC are required to be reproduced. 

2. The energy losses of hard partons in the evolving 
medium we then describe by the BDMPS mechanism [3] 
in the form presented in Ref. (23|. Assuming that par- 
tons lose energy either during the entire evolution or in 
the QGP phase only (in the hydrodynamic model), and 
varying the initial matter density profiles in the Type 2 
parametrized model, we discuss six different cases for the 
spacetime dependence of the hard parton energy losses. 
In each case, we adjust one parameter, the scale of the 
quenching power, so that the measured Raa is repro- 
duced. We discuss how Raa arises from (P(AE,E))t ab 
and which geometrical regions are probed for each case. 
In particular, we comment on the role of surface bias. 

3. Keeping all model parameters fixed from this point, we 
proceed to calculate the yield per trigger for back-to-back 
correlations for each of the 6 cases. For this purpose, we 
use a Monte Carlo simulation of the trigger setting as out- 
lined in [T3|. While each model leads (by construction) 
to almost identical Raa, we find systematic differences 
in the back-to-back yields for different associate momen- 
tum bins. We discuss both a scenario with the present 
highest trigger 8 GeV < pr < 15 GeV and the prospects 
for a higher trigger 12 GeV < px < 20 GeV. In order 
to understand what spacetime region of the evolution is 
probed in these simulations, we discuss the distributions 
of vertices in triggered events and dihadron events for 
each scenario in detail. We compare with current data 
from STAR 0, [HI and comment on the prospects of 
probing medium properties in back-to-back correlation 
measurements as well as uncertainties in the results. 



II. THE FRAMEWORK 

Our framework setup consists of three main parts: The 
hard process which is calculated in perturbative QCD 



(pQCD), supplemented by fragmentation of a hard par- 
ton outside the medium, the bulk matter evolution for 
which we either use a hydrodynamic llSfl (Type 1; 2 cases) 
or a parametrized evolution model [20( (Type 2; 4 cases) 
and the energy loss probability distribution given a hard 
parton path through the soft medium (23|. In the fol- 
lowing we describe the implementation of each of these 
ingredients in turn. 



A. The unmodified hard process 



In Ref. [2J] it has been demonstrated that leading or- 
der (LO) pQCD is rather successful in describing the 
PT-spectrum of inclusive hadron production over a wide 
range in y/s when supplemented by a v^s-dependent K- 
factor to adjust the overall normalization. This fac- 
tor parametrizes next-to-leading order effects. Since 
we are in the following only interested in ratios of Pt- 
distributions, i.e. observed yield of hadrons in A-A colli- 
sion divided by scaled yield in p-p collision, or yields per 
trigger, any factor independent of px drops out. Hence, 
in the following we use LO pQCD expressions without 
trying to adjust absolute normalization. 

The production of two hard partons k, I in LO pQCD is 
described by 



da 



AB^kl+X 



dp T dyidy 2 



Xlfi/A(xi,Q 2 )x 2 fj/B(x2,Q' 



da i: >- 



dt 



(1) 



where A and B stand for the colliding objects (protons 
or nuclei) and 2/1(2) is t ne rapidity of parton k(l). The 
distribution function of a parton type i in A at a mo- 
mentum fraction Xi and a factorization scale Q ~ px is 
fi/A( x i,Q 2 )- The distribution functions are different for 
the free protons [HI, [26[ and protons in nuclei [13, HI] . 
The fractional momenta of the colliding partons i, j are 
given by x x , 2 = ^ (exp[±yi] + exp[±y 2 ])- 

Expressions for the pQCD subprocesses ^—71 — (s, t, u) as 
a function of the parton Mandelstam variables s, t and u 
can be found e.g. in [29(. Inclusive production of a parton 
flavour / at rapidity y/ is found by integrating over cither 
yi or 2/2 and summing over appropriate combinations of 
partons, 
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da AB^f+X 

dp 2 T dy s 



d v 2 E ] 

<«>,(«) 



fcj 1 "I- O.jj 



^Xlf l / A (x 1 ,Q 2 )x 2 fj/B(x2,Q 2 ) 

0{ 



+Xlfj/A(xi 1 Q 2 )x 2 fi/B(x2,Q 2 ) 



(s, t, u)5 fk H — (s, u, t)(5/i 



dt 



dt 



dt 

da lj ^ kl 

(s, u, H -j — (I, i, u)<5/; 

at 



(2) 



where the summation (ij),(kl) runs over pairs 
<?5> 59) 59: 99i 99j 99 an( i 9 stands for any of the quark 
flavours u, d, s. 

For the production of a hadron h with mass M, trans- 
verse momentum Pt at rapidity y and transverse mass 
tot = ^/ Af 2 + from the parton /, let us introduce 
the fraction z of the parton energy carried by the hadron 
after fragmentation with z = E^/Et. Assuming collincar 
fragmentation, the hadronic variables can be written in 



terms of the partonic ones as 



rriT cosh y = zpx cosh y f and tot sinh y = Pt sinh y f . 

(3) 

Thus, the hadronic momentum spectrum arises from the 
partonic one by folding with the probability distribution 
Df—>h{ z , A 1 /) to fragment with a fraction z at a scale /i/ ~ 
Pt as 



da 



AB^h+X 



dP T dy 



E 



J 



dp T dy f 



da AB^f+X 

dp T dy f 



dzD f ^ h (z,fi 2 f )S (m T - M T (p T ,yf,z)) S(y- Y(p T ,y f ,z)) (4) 



r 



with 



and 



M T (pT,yf,z) = (zpT) 2 +M 2 t&nh 2 y f (5) 



t = \Jt 2 — z 2 , spacetime rapidity ry s = i In , cylin- 
drical radius r and in principle also on azimuthal angle 
</>, but for the time being we focus on central collisions 
only. 

The transport coefficient is related to the energy density 
of the medium as 



Pt 

Y(pt, Vfi z) = arsinh ( sinhyf 

1 tot 



(6) 



The lower cutoff z m i n = cosh y arises from the fact 
that there is a kinematical limit on the parton momen- 
tum; it cannot exceed v / s/(2coshy/) and thus for given 
hadron momentum there is a minimal z corresponding 
to fragmentation of a parton with maximal momentum. 
In the following, we use the KKP set of fragmentation 
functions Dj^h{z, fij) [30j . 



B. Energy loss of a hard parton in the soft medium 

The key quantity characterizing the energy loss induced 
by a medium with energy density e in the BDMPS formal- 
ism Q is the local transport coefficient g(r, r) s ,r) which 
characterizes the squared average momentum transfer 
from the medium to the hard parton per unit path- 
length. Since we consider a time-dependent inhomoge- 
neous medium, this quantity depends on proper time 



(7) 



with K — 1 for an ideal quark-gluon plas ma (QGP) [3l|. 
In the following, motivated by Ref. [31| we assume the 
proportionality constant K to remain unaltered in differ- 
ent phases of the medium. We treat K as an adjustible 
parameter for the following reasons: First, the energy 
loss of a parton propagating through the medium scales 
with the strong coupling a s . While we assume a s = 0.45 
throughout the energy loss calculation (and thus scale 
the q of Ref. [23[ accordingly with 0.45/0.33), the precise 
value of the parameter is not known and may be sub- 
stantially larger [Hj]. Second, we calculate energy loss 
only from the onset of thermalization but there will in 
all likelihood be some contribution from processes before 
thermalization which in principle should be accounted 
for. Third, there may also be a contribution from elas- 
tic energy loss (33[, which we do not include here, ei- 
ther. While the hard dihadron correlation pattern sug- 
gests that this is not dominant [l7j it may still contribute. 
Fourth, the temperature range probed at RHIC seems to 
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be relatively close (T < 3Tc) to the phase transition, 
thus calculating the transport coefficient in a perturba- 
tive expansion may be conceptually problematic (see [34| 
for a non-perturbative definition of q and an evaluation 
in the strong coupling limit using AdS/CFT). 
Given the local transport coefficient at each spacetime 
point, a parton's energy loss depends on the position of 
the hard vertex at ro = (xo,yo) m the transverse plane 
at r = and the angular orientation of its trajectory <f> 
(i.e. its path through the medium). To the degree to 
which the medium changes as a function of y, there is 
also a weak dependence on rapidity (in models of Type 2 
considered here). In order to determine the probability 
P(AE, E) pat h for a hard parton with energy E to lose 
the energy AE while traversing the medium on its tra- 
jectory, we make use of a scaling law [351 ] which allows to 
relate the dynamical scenario to a static equivalent one 
by calculating the following quantities averaged over the 
jet trajectory £(r) : 

/>oo 

Wc (r o ,0)= / dtm) (8) 
Jo 

and 

POO 

(g£)(r o ,0)= / m(0 (9) 
Jo 

I 



as a function of the jet production vertex ro and its angu- 
lar orientation <j>. We set q = whenever the decoupling 
temperature of the medium T = Tp is reached. In the 
presence of flow, we follow the prescription outlined in 
|36l 1371 ] and replace 



q = K-2- e 3 / 4 (p) -> K ■ 2 • e 3 / 4 (T ni " A ) (10) 

with 



T^=p(e) + [e+p(e)}-^- T (11) 



where f3± is the spatial component of the flow field or- 
thogonal to the parton trajectory. In the above two ex- 
pressions, the spacetime dependence (r] s , r, t) of pressure 
p and energy- momentum tensor T n±n± have been sup- 
pressed for clarity. 

Using the numerical results of [23|, we obtain P(AE) pat h 
for oj c and R = 2iv^/(qL) [38[ as a function of jet pro- 
duction vertex and the angle <f> corresponding to 



p(AE) path =Y.- { n / ^ 



dl(u>i) 
1 duj 



exp 



dui — 

duj 



(12) 



which makes use of the distribution toj^ of gluons emit- 
ted into the jet cone. The explicit expression of this 
quantity for the case of multiple soft scattering can be 
found in (23|. Note that the formalism of ;23J is defined 
for the limit of asymptotic parton energy, hence the prob- 
ability distribution obtained P{AE) pat h is independent 
of E. 

The initially produced hard parton spectrum and, con- 
sequently, the number of hard vertices in the (x, y) plane 
(where the z-axis is given by the beam direction) in an 
A — A collision at fixed impact parameter b, are propor- 
tional to the nuclear overlap, 



dN 



AA 



Taa(0) 



da AA^f+X 

dqlj,dy*f 



(13) 



where T^^b) is the standard nuclear overlap function 
and the cross section is from Eq. |(2]). The asterisks de- 
note the jet state before any energy losses. We define a 
normalized geometrical distribution P(xo,yo) for central 
collisions as 



[TA(ro)} 2 
T AA (0) 



(14) 



where the thickness function is given in terms of 
Woods-Saxon the nuclear density p A (r, z) as Ta(i-) = 
/ dzp A (r, z). 

Thus, we can define the averaged energy loss probability 
distribution Nl as 



Y /*27T /*00 /*00 

{P{AE)) Taa =— / d<p / dx Q / dy P(x ,yo)P(AE) path . 

27T Jo J -co J-oo 

In the following, we assume that energy loss and frag- 
mentation are cleanly separable, i.e. energy loss happens 
on the partonic level, then the hard parton emerges from 
the medium and undergoes fragmentation in vacuum. In 
practice this seems to be realized for hadrons with px > 6 
GeV (cf. [13])- 

Assuming that this condition is fulfilled, 
and that the direction (f> of an outgoing parton is not 
significantly changed, and only its energy Ei is reduced 
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by AE, we define an effective in-medium analogue for 
Eq. @, i.e. spectrum of partons which have experienced 
energy losses (not a true cross section as it is computed 
at a fixed impact parameter), 



By folding in Eqs. (|T3|) and (|T5]) . we obtain 



,~AA-^f+x 

medium 



1 



dN 



f 

AA 



dprdijf Taa(O) dprdyj 



(16) 



i~AA-^f+X 
medium 

dp T dy f 



= / dAE(P(AE)) TA/ 



d AB^f+X 

dq T dy* f d4>) ^ <%/ - y})<5(p T - (<ZT - A£))<5(0 - (17) 
I 



where again qt, J/} and 0^ are parton's kinematic quanti- 
ties before the energy loss. Inserting Eq. (fTT|) into Eq. (g]) 
yields the medium-modified spectrum 



da AA ^ h+x 

medium 



dN h AA 



dP^dy T AA (0) dP^dy 



(18) 



of hadrons originating from hard processes. 

The nuclear modification factor for central collision is 

defined as 



Raa(pt,v) 



dN\JdP T dy 



(19) 



T AA (0)daPP/dP T dy' 
With the definitions above, we now obtain it as the ratio 



,~ AA ^h+X rf pp-fc+X 

^^--s^s-'-b?*- <2o) 

as we have taken into accout the proper geometry and 
scaling already in the definition of (P(AE))t aa - 



Monte Carlo sampling of hard dihadron 
correlations 



While we are able to solve Eqs. (|4I15I17I20[) using stan- 
dard numerical multiple-dimensional integration routines 
in order to obtain R AA from the model, due to the greater 
complexity of the problem we have to rely on Monte Carlo 
(MC) simulations to obtain the yield per trigger of di- 
hadron correlations (note that this is different from [3^| 
where Monte Carlo techniques are already used to obtain 
R AA ). For sufficient statistics, the techniques should not 
lead to different results and we have verified that R AA 
can be computed also within the MC simulation with 
(within errors) identical results. 

Let us briefly explain the procedure. First, we sample 
the distribution of partons emerging from a hard ver- 
tex determined by Eq. (fTJ). This yields the parton type 
(quark or gluon) as well as the transverse momentum. 
We define randomly one of the partons as 'near side' and 



propagate it to the surface of the medium. Along the 
path, we determine lo c and (qL) by evaluating Eqs. |5J 
[9]) . The resulting values serve as input for the probability 
distribution of energy loss P(AE) pat h as determined in 

Often the plasma frequency lu c is far above the available 
jet energy and P(AE) thus extends to energies AE ^> 
£j e t- This reflects the fact that the radiative energy loss 
in [23j is derived in the limit of infinite parton energy. 
We consider a parton as absorbed by the thermal medium 
(and hence not tractable in the formalism outlined above) 
if its energy after energy loss is less than 0.5 GeV. Thus, 
in a significant number of cases the resulting outcome of 
energy loss will be absorption of the hard parton. 

We determine the actual energy loss of the near side 
parton by sampling P(AE) pat h- To find the energy of 
the leading hadron, we need the probability Pf—,h{z, /x) 
to find a leading hadron with momentum fraction z at 
scale [i. Strictly speaking, this is not the fragmenta- 
tion function as the fragmentation function yields the full 
single hadron distribution, not only the leading hadron, 
but since the trigger condition enforces on average large 
values of z, the two are virtually identical and we use 
the (normalized) fragmentation function Df—th{z,fJ,) as a 
model for P/^/ l (z,/i). Since the scale z m i n (cf. Eq. Q) 
cannot be implemented in the same way an approach 
starting from known parton properties before fragmenta- 
tion (the hadronic tut is not known at this point) we use 
a cutoff scale which is adjusted to RHIC d-Au data (see 
[i~7] ] for details). This introduces a (small) uncertainty in 
the absolute magnitude of the results once we scale the 
trigger energy upward from the measured data, as the 
cutoff scale should in principle also be altered. It does 
not alter the main result of the paper, i.e. the relative 
normalization of results for different models of medium 
evolution. 

If the resulting hadronic Ph = zpf ~ zEf fulfills the 
trigger condition we accept the event and proceed with 
the calculation of associated hadrons and the away side 
parton, otherwise we reject the event and continue the 
MC sampling by generating a new vertex. 

If an event fulfilling the trigger has been created, we de- 
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termine the kx smearing being added to the away side 
parton momentum. We sample a Gaussian distribution 
chosen such that the widening of the away side cone with- 
out a medium is reproduced. Since this is a number of 
order 1 GeV whereas partons fulfilling trigger conditions 
have frequently in excess of 15 GeV we note that this is 
a small correction. 

We treat the far side parton exactly like the near side 
parton, i.e. we evaluate Eqs. (J5J O along the path and 
find the actual energy loss from P(AE) with u> c , (qL) as 
input. If the away side parton emerges with a finite en- 
ergy, we again use the normalized fragmentation function 
Z?/^/j(z,/i) to determine the momentum of the leading 
away side hadron. If this momentum fulfills the imposed 
Pp-trigger condition for associated particle production, 
we count the event as 'punchthrough'. 

In addition, wc allow for the possibility that the frag- 
mentation of near and away side parton produces more 
than one hard hadron. We cannot simply subsume this 
in the fragmentation function as done for single hadron 
distributions as we are explicitly interested in the cor- 
relation strength between near side trigger hadron and 
other near side hadrons, thus we have to calculate sub- 
leading fragmentation processes separately. The quantity 
we need is Pf^i{z\, 2:2, /i), i.e. the conditional probability 
to find a hadron i from a parent parton / with momen- 
tum fraction z-i given that we already produced a leading 
hadron / with momentum fraction z\. In this language, 
the whole jet arises from a tower of conditional probabil- 
ities for higher order fragmentation processes. However, 
since we only probe the part of this tower resulting in 
hard hadrons, the treatment simplifies considerable. 

Moreover, since we are predominantly interested in the 
quenching properties of the medium and not in detailed 
modelling of hadron distributions inside the jet, we model 
the next-to-leading conditional fragmentation probabil- 
ity using the measured probability distribution A^zf) of 
associated hadron production in d-Au collisions [lja, [l6| 
as a function of zt where zt is the fraction of the trig- 
ger hadron momentum carried by the associated hadron. 
We include a factor 9{Ei — ^trigger — A£J — E assoc ) on the 
near side and 9(Ei — -Bpunch — AE — E assoc ) on the far 
side to make sure that energy is conserved. Note that 
associated production on the far side above the px cut is 
only possible if a punchthrough occurs. We count these 
events as 'near side associate production' and 'away side 
associated production'. 

Thus, the yield per trigger for dihadron correlations on 
the near side is determined by the sum of all 'near side 
associate production' events divided by the number of 
events fulfilling the trigger, the yield per trigger on the 
away side is given by the sum of 'punchthrough' and 
'away side associated production' divided by the num- 
ber of events (fulfilling the near side trigger condition). 
These quantities can be directly compared to experiment. 



D. Models for the medium evolution 

The medium enters the formalism through Eq. (|10[) which 
specifies the transport coefficient as a function of local 
energy density e and flow /3j_ and through the upper limit 
in the line integrals along the parton path Eqs. (J5J O 
which are terminated once the decoupling condition of 
soft matter T — Tp with Tp the decoupling temperature 
is reached. 

As outlined in Sec. HI we consider two different types 
of models for the QCD matter spacetime evolution: hy- 
drodynamics with pQCD+saturation initial conditions 
[IH and a parametrized evolution model [2(| ■ Below, we 
briefly describe the main characteristics of these models, 
more details can be found in the original publications. 
Within these models, we study the six different setups 
for the energy losses shown in Fig. [TJ 

Type I: Hydro-dynamic evolution 

In the model fl8j we use here, the hydrodynamical equa- 
tions are solved in the transverse plane assuming a lon- 
gitudinal Bjorken expansion and boost-invariance. The 
initial conditions, the initial energy densities, net-baryon 
number and formation time, are calculated from the 
pQCD + (final state) saturation model [lj|, where also 
next-to-leading order pQCD effects [4l| are effectively 
taken into account. This approach correctly predicted 
the multiplicities in central Au+Au collisions at vari- 
ous cms-energies at RHIC. To describe the multiplici- 
ties and spectra of bulk hadrons, a decoupling tempera- 
ture Tp = 150 MeV is needed [lj, HJ. The initial en- 
ergy density profile (which is correlated with the devel- 
opment of transverse flow and thus Tp) is assumed to 
be e(i"o) oc [T^(ro)] 2 . Thermalization is assumed to oc- 
cur right at formation at to = 0.17 fm/c for 5% central 
Au+Au collisions at ^snn = 200 GeV. At this early 
time, a peak energy density of 220 GeV/fm 3 is found in 
the center of the system, i.e. in a small volume. This 
value then decreases quickly as a function of time as the 
system expands. A bag equation of state is used and the 
system undergoes a 1st order phase transition from QGP 
to hadronic resonance gas at Tq — 165 MeV with a rela- 
tively long-lived mixed phase. The phase and freeze-out 
boundaries are shown by the thick lines in Fig. [1] With 
Tp = 150 MeV, the matter lifetime in the center of the 
system is about 11 fm/c and it decreases towards larger 
r. 

Setup la: Hydrodynamics. With the hydrodynamic 
evolution above, we first assume that the partonic 
jets lose energy according to Eq. (fTOf in all phases 
of QCD matter until T — Tp. In Fig. EJ top left 
panel, we show the paths of zero-rapidity parton 
jets which are produced at r = and 4.5 fm. 
Choosing a reference time of 1 fm/c, a transport 
coefficient q — 11.7 GeV 2 /fm is determined in the 
center of the medium by requiring agreement with 
the measured Raa ■ From this fit, we get K = 4.2, 
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FIG. 1: Equal temperature contour plots for the matter evolution models and energy loss setups described in the text. In the 
top left panel for the setup la, the two lines labelled T = 165 MeV show the phase boundary between the QGP and mixed 
phase, and between the mixed phase and the hadron resonance gas [l8l ]. In the setup lb, top right panel, only the QGP phase of 
the hydrodynamic evolution is considered for the partonic energy losses. In the middle and bottom panels for the parametrized 
evolution model [2(| and modifications thereof, the thick lines labelled T = 170 MeV represent the isotherm T = Tc- The 
equation of state for these models of Type II contains a cross-over transition, hence no mixed phase appears. The arrows 
indicate the path of a zero-rapidity hard parton originating from the fireball center r — and from r = 4.5 fm going radially 
outward. The values of q along the equal T contours are also indicated. 



i.e. there is some deviation from pQCD expecta- 
tions for the relation Eq. (fT0|) . 
Setup lb: Hydrodynamics+Black core. Using the 
same hydrodynamical model, we make the assump- 
tion that energy loss in the medium terminates as 
soon as T = Tc is reached, i.e. only the QGP 



induces energy loss, there is no energy loss for ei- 
ther mixed phase or hadronic gas phase. The top 
right panel of Fig. Q] illustrates this case. Since the 
initial distribution of energy density and the distri- 
bution of hard vertices both follow [Ta] 2 , and since 
the QGP exists only near the center for timescales 
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^> 1 — 2 fm/c, this has the interesting consequence 
that there is a relatively large halo of vertices sur- 
rounding the core from which a hard parton never 
encounters significant energy loss. In order to com- 
pensate for this halo and to agree with the mea- 
sured Raa the quenching power of the QGP in the 
core has to be substantial: A fit leads to K = 17.3 
(i.e. substantial deviations from pQCD expecta- 
tions) and hence at the reference time of 1 fm/c 
the transport coefficient in the medium center is 
48.75 GeV 2 /fm. Thus, this scenario is relatively 
close to geometric suppression in which there is a 
'black' region in which hard partons are always ab- 
sorbed and a 'white' region from which they always 
escape. 

Type II: Parametrized evolution 

Setup Ila: Box density. This name denotes the 
model described in [201 ] as found in a simultaneous 
fit to hadron spectra and HBT correlations. It is 
characterized by a Woods-Saxon density profile 
with a relatively small surface thickness d ws ~ 0.5 
fm, thus it somewhat resembles a box. This distri- 
bution is required by a fit to the HBT correlation 
radius R ou u as we wn l argue below in more detail, 
a sharp transition from medium to vacuum leads 
to an expanding freezeout hypersurface and this in 
turn implies peaked emission during final breakup 
of the system, i.e. a small emission duration and 
little difference between R ou t and R S ide- There is 
no microscopic justification to the use of such a 
steep profile in the initial state however. Since this 
profile is rather wide, there is no pronounced halo 
of hard vertices outside the thermalized region. 
The model gives a good description of all three 
HBT correlation radii as well as transverse mass 
spectra of pions, kaons and protons. It involves 
(non-Bjorken) accelerated longitudinal dynamics 
(for details see 20] ) . This is somewhat beneficial as 
an initial rapidity interval of ~ 4 units is mapped 
into a final interval of ~ 7 units, leading to an in- 
creased density in the initial state and hence to 
less deviation from pQCD expectations for K [3^ ]. 
The equilibration time is 0.6 fm/c. As shown in 
Fig. [TJ middle left panel, for a decoupling tempera- 
ture T F = 110 MeV the lifetime is about 20.0 fm/c. 
At the reference time of 1 fm/c, the transport coef- 
ficient in the center is 7.11 GeV 2 /fm, corresponding 
to K = 2.3. 

Setup lib: T/i-density Leaving the essential scales of 
the model [20j as above, we replace the Woods- 
Saxon density by Ta{yq) as could be expected for 
a soft matter production mechanism. With this 
modification, the model is still in agreement with 
hadron spectra and Ri ong but deviates from the 
measured R ou t and R S ide- The transverse profile is 
still wider than the distribution of production ver- 
tices, so no halo is created. With the changed den- 
sity, the evolution time is 21.5 fm/c, as is demon- 



strated by the middle left panel of Fig. [TJ At 1 
fm/c evolution time, q in the fireball center is 9.9 
GcV 2 /fm using K = 2.3. 

Setup lie: Box density+Bjorken. We keep the 
model as defined above but change longitudinal 
expansion into a boost-invariant one. This signifi- 
cantly reduces the lifetime to 17.5 fm/c, see Fig. [JJ 
bottom left panel. This is, however, still more 
than in the hydrodynamical model above which 
has a much higher decoupling temperature. The 
resulting evolution is still in fair agreement with 
R S ide, but neither R ou t nor Ri ong can be described. 
At the reference time, this leads to q = 6.85 
GeV 2 /fm in the fireball center using K = 4.0. 

Setup lid: T^-density+Bjorken. As above, but we 
change the transverse density profile into the Ta 
density. The resulting fireball lifetime is 18.5 fm/c, 
see Fig. [JJ bottom right, and at the reference time 
we find q = 9.4 GeV 2 /fm in the fireball center using 
K = 4.0. 

As discussed in [361 ] , the deviation of K from 1 is predom- 
inantly influenced by assumptions about the longitudinal 
flow and to 2nd order also determined by the magnitude 
of transverse flow. 

We note at this point that a measurement of the longitu- 
dinal expansion dynamics can be done e.g. using thermal 
photons [43j . Fixing the precise value of K would clearly 
increase the value of dihadron correlations as a tool for 
medium tomography. 

While the models presented here are far from being an 
exhaustive search through the parameter space of reason- 
able bulk matter evolution models, we believe they rep- 
resent a fair sample by including two different types of 
longitudinal dynamics, three different assumption about 
the transverse density profile and three variations in the 
decoupling parameters. 

III. SINGLE HADRON OBSERVABLES 

Let us start with a discussion of Raa and the modifica- 
tion of single hadron spectra in the model. The models 
describe Raa rather well using a single adjustable pa- 
rameter K. The quality of the description is shown for 
the hydrodynamical model in Fig. [2J it is comparable for 
the other models (not shown). 

A. High opacity saturation 

In [3!| Hi| it was argued, albeit based on simulations in 
static scenarios, that Raa for q > 5 GeV 2 /fm gradually 
loses the sensitivity to energy densities in the medium 
core, hence the observed amount of high px hadrons has 
a high probability to come from the surface. As is evident 
from the previous section, q is for some time in all of the 
models of that order (however, due to the expansion and 
subsequent dilution drops at later times). 



9 



0.8- 



• PHENIXdata 

— hydro evolution 

— hydro evolution, opacity x 1.5 

— hydro evolution, opacity x 2 




FIG. 2: Model calculation of the nuclear suppression factor 
Raa as compared to PHENIX data [IB)] for the best fit within 
the hydro model, setup la, (see text) and with additional in- 
crease of the medium opacity by 50 and 100%. In the other 
setups, the shape and magnitude of Raa are practically iden- 
tical to this figure. 



In our dynamical framework, we can test this by increas- 
ing the value of K beyond the best fit to the data. This 
directly scales the medium opacity. In Fig. [2] we show the 
best fit to Raa with the hydrodynamical models with 
parameters as given in section III Dt This gives a fair 
description of the measured pionic Raa beyond px = 5 
GeV. Increasing the opacity by 50 and 100% apparently 
leads to some kind of saturation, but the limiting curve 
for futher increases of opacity is clearly below the data 
in the region between 5 and 15 GeV transverse momen- 
tum. Thus, the conclusion is that if the full dynamics 
is taken into account, saturation of Raa with respect to 
increasing q and dominance of surface emission is not yet 
reached. For similar conclusions in a different framework 
see also I46II. 



B. The geometry underlying Raa 

We can gain further insight into the question of surface 
vs. volume dominated emission by studying the origin of 
trigger hadrons from the MC simulations. These contain 
the same information as Raa, albeit in more differential 
form at a given scale (i.e. the trigger momentum). This 
is shown in Fig. [3] 

As evident from the figure, while there is some degree of 
surface bias and little emission is found from the fireball 
core, the spatial region probed by a single hadron distri- 
bution extends deep into the medium. The only scenario 
where this is not quite the case is the 'black core' hy- 
drodynamics. Here, a clear trend to surface emission is 
visible, although still the core is not completely black. 
It is evident that the degree of surface bias is model de- 
pendent. We find consistently that in scenarios in which 
the initial distribution of matter is rather wide (and cor- 



respondingly the local q can be smaller to achieve the 
same Raa) the degree of surface bias is reduced. There 
is also some sensitivity to the underlying density profile 
visible in the comparison of the box and Ta density sce- 
narios. This gives some confidence that once the away 
side parton is considered, its sensitivity to (P(AE))T r 
will be able to distinguish two scenarios which which are 
characterized by the same {P(AE))t ab but have differ- 
ent distributions of matter. 



C. The role of transverse flow 

Let us at this point briefly remark on the role of trans- 
verse flow. In [17| it was pointed out that this is a cru- 
cial effect for long pathlength and increases the trans- 
parency of the medium. While the transverse geometry 
does not change significantly for timescales of 2-3 fm af- 
ter equilibration since transverse flow takes some time 
to develop, the change of transverse geometry by trans- 
verse flow clearly is an issue for longer timescales (which, 
according to Fig. [3] are frequently probed even for sin- 
gle hadron observables when the vertex lies close to the 
fireball center). 

It is sometimes argued that the effect of transverse flow 
cancels for radiative energy loss with a quadratic depen- 
dence on pathlength. This argument goes as follows: 
Suppose we have a homogeneous medium with expands 
with the velocity vt in radial direction. The medium 
density (and hence the transport coefficient) drops like 
1/(vtt) 2 due to the transverse expansion. However, the 
freeze-out surface also moves outward with vt, and hence 
there is an additional pathlength ~ vtt the particle has 
to go through the medium. For quadratic pathlength de- 
pendence of energy loss, this implies additional energy 
loss ~ (vtt) 2 which would just cancel the effect of the 
dropping transport coefficient in Eq. ([5]). 
What is missing in the argument is that an actual 
medium is not homogeneous, and therefore the freeze- 
out hypersurface does in general not expand with the 
flow lines. Only in the case of a homogeneous box den- 
sity profile is this the case, for any realistic density profile 
the freezeout surface expands much more weakly or even 
shrinks in time. This can clearly be seen e.g. in Fig. 2 
of Ref. [18|. However, especially in scenarios where the 
freezeout radius shrinks with time, transverse flow has 
an enormous impact — not only is the density dropping 
with (vtt) 2 but there is also a systematic shortening of 
the in-medium paths. Thus, transverse flow can, con- 
trary to the naive expectation, have a significant impact 
on the medium transparency. 



D. Average quenching properties 

In Fig. 3] we show the average energy loss per unit time 
(unit pathlength) as a function of time for a quark start- 
ing from the center r = of the fireball and from the 
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FIG. 3: Probability density of finding a parton production vertex at (a;, y) given a triggered event with 8 GeV < pr < 15 
GeV for different spacetime evolution scenarios (see text). In all cases the near side (triggered ) hadron propagates to the —x 
direction, hence the y — (—y) symmetrization. Countours are at linear intervals. 
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typical emission region (the maximum seen in Fig. [3l 
r w 4.5 fm) propagating radially outward. We obtain 
this quantity via 

/>oo 

(AE) = / P{AE) T AEdAE (21) 
Jo 

with P(AE) T obtained with the help of the line inte- 
grals Eq. ([8]) and Eq. ([9]) with the upper integration limit 
changed into r. Interestingly enough, the total average 
energy loss in all scenarios is between 20 and 23.5 GeV 
for a quark (between 34 and 37 GeV for a gluon), al- 
though evolution of the loss per unit pathlength is quite 
different in the different models. One may speculate that 
this is caused by the fact that all scenarios describe Raa 
by construction. Since the measured value of Raa can 
be understood by the observation that about 80% of all 
partons are not observed in the perturbative region and 
since most partons originate from the region around the 
medium core, the similarity of the average energy lost 
presumably also ensures that the probability distribu- 
tion of energy loss for the typical parton from the cen- 
tral region is comparable and hence a similar fraction of 
absorbed partons from this dominant central region ob- 
served. However, the origin of the observed minority of 
partons is quite different in all cases. 
Qualitatively, all curves exhibit the same shape. First, 
there is a strong rise: As pathlength increases, it allows 
for decoherence of softer and softer quanta (parametri- 
cally the decoherence length of a radiated quantum with 
momentum qx transverse to the hard parton and energy 
u> goes like Td ec ~ iolk %\ h owever soft quanta are more 
likely to be radiated) [23]. In a static, homogeneous 
medium this feature leads to the quadratic pathlength 
dependence of radiative energy loss. In a dynamic evo- 
lution, there is eventually a turnover as the density of 
the medium is decreased, either because volume expan- 
sion dilutes the medium over time or because the parton 
reaches the thinner outer layers of the transverse density 
profile. This turnover point is model dependent: for the 
hard partons produced in the fireball center (left panel), 
the average energy loss per time peaks at 1-3 fm for the 
hydrodynamic scenarios Ia,b, and at 2-4 fm for the se- 
tups Ila-d. This implies that if a parton can escape the 
medium during the first 1-2 fm/c evolution time it will 
never undergo substantial energy loss (this is not quite 
true in the 'black core' scenario though). Thus, the halo 
region where partons are not strongly medium modified 
may be expected to reach into the medium at least 1-2 
fm beyond the position of the freezeout hypersurface; de- 
pendent on the density profile at the outer edge of the 
medium even further. 

The later flattening of the curve, most pronounced in 
the hydrodynamical evolution after about 4 fm/c is re- 
lated to the slowdown of the expansion rate in the mixed 
phase when the pressure of the medium vanishes and ac- 
celerated transverse expansion turns to constant trans- 
verse expansion. This feature is less pronounced in 



the parametrized evolution scenarios as these employ a 
crossover transition where a soft point in the EOS is 
reached but the transverse acceleration never drops to 
zero. 

In the 'black core' scenario, the average energy loss 
of partons released from the center gets as high as 18 
GeV/fm and even in the least dramatic box density case 
it reaches up to 4 GeV/fm. Thus, it is safe to conclude 
that the medium is on average extremely black — prop- 
agation of partons from the medium core to the surface 
would take of the order of 25 GeV parton energy, thus 
the typical parton energy before energy loss and frag- 
mentation would have to be ~ 40 GeV to form a 8 GeV 
hadron, a scale far above typical parton energies available 
in significant numbers. We therefore conclude that fluc- 
tuations around the average energy loss must be large. In 
a situation in which partons are absorbed by the medium 
on average, fluctuations will open up the possibility of the 
parton being able to penetrate the medium and hence in- 
creases the transparency somewhat. Let us explore this 
by studying a probabilistic representation of the quench- 
ing process in the model. 



E. Averaged energy loss probabilities 

We show the geometry-averaged energy loss probabilities 
for quarks (P(AE))t aa (see Eq. ([15))) for the different 
scenarios in Fig. [5j This quantity reinforces our conclu- 
sion from the average energy loss: Typically quenching 
is substantial, but there are strong fluctuations. The 
probability distributions exhibit long tails extending out 
above 100 GeV energy loss, but there is also a large es- 
cape probability of 0.26 for the hydrodynamics case, 0.3 
for the black core case and 0.24 for all other scenarios. 
Note that all parametrized evolutions lead to virtually 
the same averaged energy loss probability — they would 
be indistinguishable even by a 7-hadron correlation mea- 
surement as outlined in [8j. 

Taking gluons into account, there is little actual energy 
loss observed in the model: About 15-20% of the par- 
tons escape without energy loss (either because they are 
created outside the medium or due to fluctuations in the 
energy loss probability) , only about 5-8% of partons con- 
tribute to the hard hadron spectrum after undergoing 
some energy loss and 75-80% of partons are absorbed 
in the medium and thermalize. Thus, the information 
about the medium is predominantly carried in the ratio 
of absorption to transmission, not in the average energy 
loss of observed hadrons. 

This has been observed in [Tt], [39[ already and seems 
to be characteristic for the BDMPS energy loss in the 
formalism of [23J] in the RHIC energy range as long as 
a realistic distribution of pathlcngths is taken into ac- 
count. Thus, if BDMPS energy loss is realized in nature, 
one has to go even beyond 7-hadron correlations (which 
reflect the [T^(r)] 2 profile of the hard vertices) to gain 
sensitivity to details of the medium density distribution. 
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FIG. 4: Left panel: Average energy loss per time for a hard quark released in the fireball center r = propagating towards 
the surface for the different setups studied. Right panel: Same quantity but for a quark propagating radially outward from the 
maximum of emissivity seen in Fig. [3] i.e. about 4.5 fm from the center. Notice the difference in the horizontal scales. 
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FIG. 5: Geometry-averaged energy loss probability 
{P(AE))t aa of quarks for different spacetime evolution sce- 
narios (T_A-density, box density Bjorken and Ta density 
Bjorken are not shown, they are indistinguishable from the 
box density scenario). 



In the following, we explore the possibility of doing this 
in hard back-to-back correlations of hadrons. 



ditional suppression of the away side is about of the same 
order as the near side suppression — roughly 4 in 5 of 
partons back to back with a valid trigger are quenched. 
This is an interesting observation in itself, as it clearly 
demonstrates that the systematic difference in path- 
length between near and away side is important. This 
essentially rules out a purely geometric interpretation of 
jet quenching where partons born in a 'black' region are 
always absorbed and partons born in a 'white' region al- 
ways survive — in such a scenario there would be no addi- 
tional absorption of away side partons and the yield per 
trigger (modulo fragmentation) would be of order one. 
But it also places strong constraints on the time cutoff of 
quenching — if the medium becomes transparent due to 
the volume expansion too soon, there is no time for the 
away side partons to pick up additional energy loss due 
to their longer in-medium path. Unfortunately, in order 
to make these statements quantitative, a microscopic de- 
scription of jet energy loss including the full geometry 
(vertex distribution in Fig. [3]) is required. Since we have 
limited ourselves in the present investigation to BDMPS 
radiative energy loss, we will not explore this interesting 
possibility further here but rather leave this to a subse- 
quent publication. 



IV. DIHADRON OBSERVABLES 

Let us at this point remark that back-to-back dihadrons 
are a rare event and do not reveal the typical situation 
of a parton pair emerging from a hard vertex but rather 
a highly unlikely coincidence. About 4 out of 5 partons 
potentially leading to a hard trigger above 8 GeV are ab- 
sorbed by the medium, and the yield per trigger on the 
away side is on the order of 2% [lU, [l6| , thus there is 
massive additional suppression. However, much of this 
is due to the low probability of a hard away side parton 
to fragment into a hard hadron. Taking this effect into 
consideration by comparing with the d-Au data, the ad- 



A. The geometry of dihadron correlations 

We show the probability density of vertices in the (x, y) 
plane leading to a near side trigger hadron above 8 GeV 
and an associate away side hadron with 4 GeV < px < 6 
GeV in Fig. [51 It is immediately obvious that the distri- 
bution is very different from the distribution of vertices 
for single hadron observation shown in Fig. [3J First, the 
dihadron distributions are much wider in ±y direction, 
indicating the importance of the periphery where both 
near and away side parton have a short in-medium path 
(or the halo where the production vertex lies outside the 
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FIG. 6: Probability density for finding a vertex at (a;, y) leading to a triggered event with 8 GeV < pr < 15 GeV and in 
addition an away side hadron with 4 GeV < pr < 6 GeV for different spacetime evolution scenarios (see text). In all cases the 
near side hadron propagates to the —x direction. 
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medium) . This is most clearly seen for the black core sce- 
nario where there is almost complete repulsion of events 
from the dense core to the periphery. 
The second observation is that while the single hadron 
vertex distributions are typically centered around x ~ 
—4... — 5 fm, the dihadron distributions roughly center 
around x ~ where both near and away side have similar 
pathlengths. Clearly there is also pronounced sensitivity 
to the medium density — in both Ta density scenarios 
the distribution is repelled from the center of the fireball 
by the peak in medium density — no such trend is seen 
for the box densities (where no pronounced peak in the 
fireball center is present). Thus, just based on the ob- 
served geometry one would conclude that there is signifi- 
cantly more sensitivity to medium properties in dihadron 
correlations than in single hadron suppression. 



B. Comparison with STAR data 

We show the results of the model calculation for near and 
away side yields per trigger for a trigger of 8 GeV < px < 
15 GeV as a function of associate hadron momentum bin 
in comparison with the data (TBI . [l6j for central collisions 
in Fig. [3 

Within errors, the near side yield per trigger is described 
by all the models well. There is no significant disagree- 
ment among the models. This is not very surprising — 
as we have seen above and remarked in [T3|, about 80% 
of all near side partons emerge from the medium without 
having experienced energy loss. Thus, it is not expected 
that energy loss is able to modify the next-to-leading 
fragmentation of the trigger parton significantly. 
The model calculations appear significantly more differ- 
ent if we consider the away side yield. Here, results for 
the 4-6 GeV momentum bin differ by almost a factor 
two. However, none of the model calculations describes 
the data in this bin. This is in fact not at all surprising as 
below 5 GeV the inclusive single hadron transverse mo- 
mentum spectra are not dominated by pQCD fragmen- 
tation and energy losses but, rather, by hydrodynamics 
possibly supplemented with recombination pTl |48| type 
phenomena, see Fig. 9 of Ref. 18]. For this reason, 
the ratio Raa at px < 5 GeV cannot be expected to 
be described by pQCD fragmentation and energy losses, 
either. However, the yield of hadrons associated with a 
given trigger must reflect the structure of the underly- 
ing event. Therefore, the (uncorrelated) recombination 
of thermal partons cannot be responsible for the dis- 
crepancy. Rather, in the language of [48|, thermal + 
shower recombination processes are likely candidates for 
the missing contribution, likewise a possible distortion of 
the underlying hydrodynamical flow by the thcrmaliza- 
tion of lost energy. Both these contributions are expected 
to be small in the 6+ GeV momentum region. Given that 
our model at present incorporates only hadron produc- 
tion by fragmentation, it would be a mere coincidence if 
a good agreement of the low-energy bin on the away side 



were obtained between the data and our calculation. We 
conclude that our model cannot offer a reliable prediction 
of away side hadronic yields in this bin. 
This is clearly unfortunate, as the model results are con- 
siderable closer to the experimental result in the 6+ mo- 
mentum bin on the away side and hence our ability to 
discriminate between different models is reduced. Since 
at this large transverse momenta the pQCD fragmenta- 
tion + energy losses dominate the singe hadron spectrum 
(again see Fig. 9 of Ref. [Hj]), we expect that the model 
is able to give a valid description of the relevant physics 
in this bin: Not only is Raa well described by the data, 
but also the contribution of recombination processes to 
the yield is expected to be small [47j . Thus, as it stands, 
only the black core scenario can be ruled out by the data, 
the box density with Bjorken expansion seems strongly 
disfavoured but still marginally acceptable. In a sense 
this is certainly reassuring, as this indentifies the sce- 
nario least likely to be realized in heavy-ion collisions — 
the black core scenario exhibits strong deviations from 
pQCD expectations for the energy loss and there is no 
a priori reason that the mixed phase or the hadron gas 
should not contribute to energy loss. 
Nevertheless, while there are indications that the other 
scenarios show sizable differences in the momentum spec- 
trum of away side hadrons, the present data is not suf- 
ficient to make a distinction. There are in principle two 
ways to overcome this problem. At the price of intro- 
ducing additional model dependence, one might include 
recombination processes into the simulation. In this way, 
comparison to the 4-6 GeV momentum bin would be pos- 
sible. Alternatively, one can address the question if more 
leverage in px would improve the question and hence if 
improved experimental conditions will allow tomography. 
We have chosen to follow the latter path. Towards this 
end, we study in the following a situation where the trig- 
ger momentum is increased to 12 GeV < px < 20 GeV 
and consequently more momentum bins in the pQCD re- 
gion become accessible. 



C. Partonic and hadronic momentum spectra for 
pr > 12 GeV trigger conditions 

We present the distribution of away side parton momenta 
given a hard triggered hadron on the near side in Fig. [5] 
before (left panel) and after (right panel) energy loss due 
to the medium. The spread of the distribution before en- 
ergy loss is a measure for the amount of energy loss (as 
compared to transmission or absorption) induced by the 
scenarios; or equivalently, the similarity of the curves re- 
states the fact that the geometry-averaged distributions 
(P(AE))x AA are rather similar. 

Considering the distribution after energy loss (and dis- 
regarding again the dominant absorption contribution), 
there is a sizable shift of the spectral distribution towards 
lower momenta, but this shift is different for the different 
models. Clearly, the black core model is most extreme, 
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FIG. 7: Yield per trigger on the near side (left panel) and away side (right panel) of hadrons in the 4-6 GeV and 6+ GeV 
momentum bin associated with a tr igg er in the range 8 GeV < pr < 15 GeV for the different models of spacetime evolution as 
compared with the STAR data [15l, Il6|. The individual data points have been spread artificially along the x axis for clarity of 
presentation. 



0.2 



0.15 



0.1 



0.05 



hydrodynamics 
black core 
box density 
T A density 
box density, Bjorken 
T A density, Bjorken 



0.2 



0.15 



-Oh al- 



io 20 

P T [GeV] 



30 



40 



HI 



0.05 



» hydrodynamics 

• black core 

• box density 
. T A density 

• box density, Bjorken 
T A density, Bjorken 



10 



20 
P T [GeV] 



30 



40 



FIG. 8: Conditional probability distribution P(pr) to find momentum px for the away side parton given a triggered near side 
hadron in the range 12 GeV < Pt < 20 GeV before (left panel) and after (right panel) away side energy loss due to passage 
through the medium for different spacetime evolutions. Here, we only consider away side partons not absorbed by the medium. 



but differences of the order of 20-30% are also seen in the 
other models. 

The distribution after fragmentation into hadrons in bins 
of 2 GeV width in the perturbative region is shown in 
Fig. [9] for the near side (left panel) and away side (right 
panel) . It is again apparent that within errors all models 
agree in the expected near side yield. The momentum 
spectrum of the away side exhibits considerably more 
structure. Several of the scenarios can now be clearly 
told apart in bins in the perturbative region. For ex- 
ample the Ta and the box density which have virtually 
identical (P(AE))t aa (cf. Fig. [5]) show almost a factor 
two difference in the 10-12 GeV momentum bin; Ta and 
Ta Bjorken can be told apart in the 6-8 GeV momentum 
bin (provided enough experimental precision is achieved) . 

It is evident from the analysis that having a larger lever- 
arm in momentum is clearly beneficial — as apparent 
from the figure, the momentum distribution of away side 
hadrons after energy loss with (P(AE)) trigger is char- 



acteristic of the scenario, and although the differences 
induced by the different geometry and expansion pattern 
are not factors of 10, they may reach as much as 50%. 



D. Changes in geometry with trigger and associate 
energy 

In Fig. [TU] we investigate to what degree the region in 
the transverse plane probed by the dihadron correlation 
is changed with increased trigger momentum and/or as- 
sociate cut. We do this at the example of the hydrody- 
namical evolution. 

In the left panel, we show the distribution of vertices 
leading to an away side hadron between 4 and 6 GeV, but 
for a trigger hadron above 12 GeV. Comparing Figs. ITD1 
andO there is (given the limited statistics) no significant 
difference between the figures. 

However, going to higher associate hadron momentum 
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FIG. 9: Yield per trigger on the near side (left panel) and away side (right panel) of hadrons in the 4-6 GeV, 6-8 GeV, 8-10 GeV, 
10-12 GeV and 12+ GeV momentum bin associated with a trigger in the range 12 GeV < pr < 20 GeV for the different models 
of spacetime evolution. The individual data points have been spread artificially along the x axis for clarity of presentation. 
Note that the last bin extends from 12 GeV up to the pr of the trigger hadron and is thus considerably wider than the previous 
bin, explaining the upward turn of some spectra. 
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FIG. 10: Probability density for finding a vertex at (x, y) leading to a triggered event with 12 GeV < Pt < 20 GeV and 
in addition an away side hadron with 4 GeV < Pt < 6 GeV (left panel) and 8 GeV < Pt < 10 GeV (right panel) for the 
hydrodynamical evolution model. In all cases the near side hadron propagates to the —x direction. 

between 8 and 10 GeV, we find that the distribution is mentum scaling, 
somewhat more pushed out of the fireball center. Appar- 
ently, in this case even small energy loss is disfavoured. 
But all in all, the influence of the spacetime distribution 
of matter on the probability distribution of vertices is 
greater than the influence of trigger and associate mo- 
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FIG. 11: Model calculation of the away side yield per trig- 
ger for an 8 < pr < 15 GeV hadron trigger in two different 
momentum bins as compared to STAR data [TEL [l(| for the 
hydro model, setup la, (see text) and with additional increase 
of the medium opacity by 50 and 100%. 

E. High opacity saturation? 

Since hard dihadron correlations select a special class of 
events with an underlying vertex distribution very dif- 
ferent to the distribution underlying any single hadron 
observable, it is useful to investigate the question if any 
saturation of the yield with the quenching power of the 
medium is reached. Note that there is no a priori rea- 
son to expect the density at which single and dihadron 
observables saturate to be the same, as dihadron observ- 
ables probe the quenching of rare fluctuations in the en- 
ergy loss probability distribution. 

As in the case of Raa, we can test this by increasing the 
value of K beyond the value obtained from a fit to Raa- 
We show the resulting reduction in the away side yield 
per trigger in Fig. [TT] calculated within the hydro model 
setup la. There is some flattening of the curves observed 
but no saturation of the yield within the opacity given 
by the model. This observation is consistent with the 
information from Fig. [5]where it is clearly seen that apart 
from the black core scenario hard dihadrons originate 
from the medium center and are not yet pushed out to 
the periphery by unpenetrable central densities of the 
medium. The fact that saturation of Raa appears to be 
reached earlier than saturation of the away side yield in 
dihadron correlations, in particular in the 6+ momentum 
bin, is also consistent with the observation that single 
hadron observables show generally more surface bias (see 
Fig.©. 

V. DISCUSSION 

We have studied the possibility of doing jet tomography, 
i.e. discriminating structures of the density distribution 



and the expansion pattern of the medium created in ul- 
trarelativistic heavy-ion collisions using back-to-back cor- 
relations of hard hadrons. 

We found that due to the different geometry enter- 
ing the averaged energy loss probability distributions 
(P(AE))t ab and (P(AE))t t there is non-trivial infor- 
mation in the dihadron yield per trigger beyond what 
is constrained by Raa- We have explicitly shown that 
even the current data are able to rule out a somewhat 
more extreme scenario of jet quenching and that a greater 
lever-arm in away side hadron momentum offers the pos- 
sibility to discriminate scenarios which lead to the same 
(P(AE))t ab and are hence even in principle indistin- 
guishable by either Raa or 7-hadron correlation mea- 
surements. 

The requirement of having a large number of momen- 
tum bins to get better discrimination however makes this 
method in all likelyhood more suitable for the LHC where 
hadrons up to 100 GeV momentum are expected to be 
observed regularly and thus a large number of momen- 
tum bins could be measured with great precision far in 
the spectral region where pQCD and vacuum fragmenta- 
tion is applicable. 

Given the large possible parameter space of medium evo- 
lution models, it is unlikely that any form of jet tomog- 
raphy alone will yield a complete characterization of the 
medium density. It seems rather that a multi-pronged ap- 
proach, i.e. simultaneously measuring 7-hadron correla- 
tions (and hence (P(AE))t ab for quarks), back-to-back 
correlations (i.e. (P(AE))Tr), thermal photons (sensi- 
tive to the longitudinal expansion [43|) and a reaction 
plane analysis of Raa [S] introducing a handle on the 
systematic variation of the in-medium pathlength dis- 
tribution will be a suitable tool to extract tomographic 
information using hard probes. 
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